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The trapped w-modes of stars with a first order phase transition (a density discontinuity) are 
computed and the excitation of some of the modes of these stars by a perturbing shell is investigated. 
Attention is restricted to odd parity ("axial") m-modes. With R the radius of the star, M its mass, 
Ri the radius of the inner core and Mi the mass of such core, it is shown that stars with R/M > 5 
can have several trapped w-modes, as long as Ri/Mi < 2.6. Excitation of the least damped ii;-mode 
is confirmed for a few models. All of these stars can only exist however, for values of the ratio 
between the densities of the two phases, greater than ~ 46. We also show that stars with a phase 
transition and a given value of R/M can have far more trapped modes than a homogeneous single 
'j ' density star with the same value of R/M, provided both R/M and Ri/Mi are smaller than 3. If the 

, phase transition is very fast, most of the stars with trapped modes are unstable to radial oscillations. 

■ We compute the time of instability, and find it to be comparable to the damping of the w-mode 

^S) ' excited in most cases where w-mode excitation is likely. If on the other hand the phase transition 

J, , is slow, all the stars are stable to radial oscillations. 

\o ■ 

; I. INTRODUCTION AND OVERVIEW 

. The w-modes of a fluid star ||l|-^ represent pulsations of the spacetime geometry of the star, during which the fluid 
' remains almost at rest. In general these modes damp much faster than the fluid modes of the star (although very 
', slowly damped w-modes do exist for very compact stars Q). This raises the possibility that, in case the w-modes are 
f/^ ■ astrophysically excited, they might carry an important part of the gravitational energy radiated by the pulsating star. 
^ Although stars of any compactness seem to have w-modes , the damping of these modes is so fast for Newtonian 
' stars that even if one of the modes is excited, the star only has time to oscillate a tiny fraction of one period before 
' the vibration dies out. The w-modes are thus largely irrelevant compared to the fluid modes for noncompact stars. It 
O is only when the star becomes sufficiently compact that the least damped w-modes have damping rates slow enough 
that the star might undertake a few complete periods of oscillation at the w-mode frequency. 
1 Recent investigations [p|-|l4| have started to address the conditions under which w-modes of compact stars are 
excited. Nearly all of them have shown that only unrealistically compact stars, with a radius to mass ratio R/M < 3, 
^ ' have their least-damped w-modes strongly excited. For neutron stars (or any other known stars), which have R/M > 5, 
the w- modes are at best weakly excited. The star pulsates essentially through its fluid modes. It is important however 
to point out some limitations of these studies. On the one hand they consider pulsations of static, spherically symmetric 
^ . stars, i.e. they neglect effects due to the rotation of the star. On the other hand, they consider sources of perturbation 
■ " ' that probably do not mimic entirely the main aspects of a astrophysically realistic dynamical situation that might 
leave behind pulsating remnants: the collision of two neutron stars, a supernova explosion or the gravitational collapse 
of a massive star. There is also the issue of the equation of state used to model the matter of the star. The above 
studies restricted their attention almost exclusively to polytropic equations of state or equations of state in which the 
density is a smooth function of the radial position of the matter. While this is reasonable for typical neutron stars, 
the w-modes depend sometimes dramatically on the equation of state and their excitation will likewise be affected. 
In particular the effect of a flrst order phase transition (a discontinuity in the density) appears not to have been 
considered before in connection with w-mode studies. 

In this work, we use a rather different, even if somewhat unrealistic equation of state, to show that some stars with 
R/M > 5 have a w-mode spectrum very different from the spectrum of a star described by the usual equations of 
state. We then go on to show that the least damped of these w-modes are excited by matter moving on a spherical 
shell surrounding the star. As in previous works |^,|o| we consider the unperturbed star to be static and we regard 
the shell as a perturbation of the spherically symmetric spacetime. Because the main intention of this work is to show 
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the effect of the choice of the equation of state on the w-modes and on their excitation, we restrict attention to odd 
parity (axial) perturbations and hence to the odd parity quasi-normal modes of the star, which are all w-modes. 

We note that a partial reason why the w-modes of stars with R/M > 5 are barely visible and hence hardly relevant, 
is because they are very quickly damped. The star simply does not have time to oscillate. On the contrary stars 
with R/M < 3 have w-modes that are either mildly or, for some equations of state, slowly damped. Usually only for 
these stars are the w-modes significantly excited, carrying sometimes the same or more energy than the excited fluid 
modes. 

The class of slowly damped modes (sometimes called trapped modes) was first discovered by Chandrasekhar and 
Ferrari Q. They looked at odd parity oscillations of spherically symmetric stars, which are governed by a Schrodinger- 
like equation. Drawing from studies in atomic physics, they pointed out that a necessary condition for existence of a 
trapped w-mode is that the potential of the Schrodinger equation has a minimum followed by a maximum (a slightly 
more restrictive necessary condition has been proposed in pst). This is not, however, a sufficient condition: only 
when the depth and width of the potential well is high or large enough will one (or more) trapped w-mode(s) exist. 
Chandrasekhar and Ferrari then went on to determine the trapped ui-modes of homogeneous stars and found that 
they exist only for stars with 2.25 < R/M < 2.6. The value R/M — 2.25 is the highest compactness such stars can 
reach. 

Early on it was believed that in order to have trapped w-modes, a star should have values of R/M smaller than 3. 
But recent work by Rosquist jl^] and by Karlovini, Rosquist and Samuelsson ||l7|-[T9[| suggests that there are equations 
of state for which stellar models with R/M > 3 having trapped w- modes exist. They studied odd parity oscillations 
of spherically symmetric stars and showed that certain equations of state lead to a potential (in the Schrodinger-like 
equation governing these oscillations) having one or more wells, even when the star has R/M > 3. If one or more of 
the wells are deep and broad enough, slowly damped w-modes should exist. And if they exist, it is possible that the 
star might oscillate in one or more of these modes when perturbed by a plausible astrophysical source. 

In this work we consider the simplest of the equations of state studied in p^-pj[|. The static star has an inner 
core of constant density p+, surrounded by an envelope of constant density p- such that p- < p+. We will call 
these double-density stars (DDS). As long as we allow the jump in the density to be sufficiently high, the following 
argument of Rosquist jl^ suggests strongly that some DDS models with arbitrarily high values of R/M might have 
slowly damped w-modes. Let us consider a very compact star of constant density p+, with radius Ri and mass Ali 
such that Ri/Mi < 2.6. This star has trapped w-modes. Let us then transform a thin layer of width 6, located at the 
surface of the star into an envelope with some different but also constant density p_ < p+, in such a way that the 
mass of the whole star is kept approximately constant. The radius of this new star is i? « i?i(l -I- 2ip+/ pS / RiY^^ . 
It is then clear that choosing p+/ P- » Ri/5 » 1 will lead to R/M » Ri/Mi, since we kept M w Mi. But we 
expect the new DDS star to have trapped modes similar to the old star since all we did was add an atmosphere of 
negligible density to the original star. Additional evidence that this argument should be correct is given in |p7|| . 

It is important to ask what is the physical plausibility of DDS stars (and ultimately of the stars considered in 
[p^-p^). A realistic neutron star can be seen as being formed by several layers. The density in these layers increases 
(usually but not always) continuously but over a very short length scale, as we approach the center of the star. It 
increases from p < lO^g/cm^, in the surface region to p > lO^^g/cnv^ in the inner core. A DDS star could then be 
viewed as modeling one of the abrupt density jumps in a very crude way: by supposing it to be discontinuous. A 
second case in which a DDS equation of state might be appropriate is during the gravitational collapse of a massive 
star. There it might be possible that some of the layers experience a jump in the density, eventually well represented 
by a DDS star. But this dynamical situation would likely be of little relevance to the issue of w-mode excitation 
since the time of collapse would be much smaller than the period of one of the least damped w-modes. A third, 
much more realistic hypothesis, is to regard a DDS as a model of a star with a first order phase transition. Some 
of these transitions do occur in realistic neutron stars. Directly below the surface of the neutron star, lies the outer 
crust, a solid region of heavy nuclei in equilibrium with a relativistic degenerate electron gas where the density ranges 
from lO^g/cm^ to 4 x lO^^g/crrv^. At densities p > W'^g/cm'^, the electrons combine with bound nuclear protons to 
form nuclei with a higher number of neutrons through inverse (3 decay. This neutronization leading from one nuclei 
to another is a first order phase transition. At p « 4 x 10^^ g/cm'^ neutron drip sets in. The region bellow the 
outer crust, the inner crust, with densities in the range 4 x lO^^g/cm^ < p < 2 x lO^^g/cm^ is thus a mixture of 
neutron-rich nuclei together with a superfluid neutron and electron gas. The transition from the solid outer crust to 
the liquid inner crust can also be seen as a first order phase transition. In both of these phase transitions however, 
the density jump 7 — p-^-/p^ is quite modest, w 3% at the most. Very massive stars might have regions where 
p > 10^^ g/cm^. In those phase transitions from a phase with nucleonic matter to quark matter may occur. Other 
exotic phase transitions are also predicted at these very high densities: kaon, pion and hyperon condensation. These 
more speculative phase transitions, might lead to jumps 7 of up to 10. For a review about phases in neutron stars 
and corresponding bibliography from the point of view of nuclear physics see ||2^ . 

Two-phase stars have been modeled both with Newtonian and general relativistic gravity, starting with the pi- 
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oneering work of Ramsey pi| in 1950. He obtained the equilibrium configurations of Newtonian stars with two 
incompressible layers of constant density and showed that some of these stars might be dynamically unstable to radial 
oscillations. Ovakimova considered the general relativistic analog of Ramsey's star and derived the corresponding 
stability criteria. Stability criteria for stars whose phases are not of constant density are scarce, except in the case 
when the core of the denser phase is very small. In this case, Lighthill generalized Ramsey's stability criteria 
to two-phase Newtonian stars with phases described by any equation of state, and showed that all such stars are 
unstable if the jump in the phase density is bigger than 3/2 and are stable otherwise. Seidov extended Lighthill's 
results to general relativistic two-phase stars and showed that, if the core of the denser phase is very small, general 
relativity increases the range of stability of a two-phase star (see also [^). Reviews and recollections by a pioneer in 
the field of two-phase stars are |^,^. See also p^-|30| for several models of equilibrium two-phase stars. 

In this work we will look at a DDS as representing a star with a first order phase transition. We will also make 
the additional assumption that each of the layers of the star is incompressible. Since the main point of this work is 
to show that stars with R/M > 5 have slowly damped ui-modes that can be excited, it is important to distinguish 
between two limiting types of phase transition (a similar distinction is made in ||3l|]). A slow phase transition is one 
that happens on a time scale much longer than any of the least damped w-modes. As far as tu-mode excitation, these 
stars behave as if there were no phase transition at all and the star where made of two different substances in contact 
along a common boundary. A fast phase transition is one that happens on a timescale much shorter than the period 
of the least damped w-modes. 

Each type of star will behave the same way under odd parity perturbations but differently under even parity 
perturbations. This is manifest for example if we consider radial oscillations. If the phase transition is very slow then 
radial oscillations can not exist, the star being incompressible. If on the contrary the transition is fast then the star 
can oscillate radially |^,^,^,^ . During this oscillation some of the matter in one of the phases will be converted 
into the other phase and vice- versa. It would appear that which type of star we consider, is not relevant if our interest 
is only in odd parity perturbations. However stars with a fast phase transition might in some cases become unstable 
to the (even parity) radial oscillation, as we mentioned above. If this instability time is much longer than damping 
time of an excited odd parity w-mode, then which type of star we consider is indeed irrelevant. If on the other hand 
the instability time is comparable (or much smaller) than the damping time of the w-uiode it is important to make 
explicit which type of star we are considering. In this later case, it does not make sense to talk of excitation of 
w-modes of a two-phase DDS, while it makes sense to investigate this excitation if the phase transition is slow, for 
then the star will not be unstable. 

We show that stars with R/M > 5 have trapped modes as long as the compactness of their core, as measured by 
the ratio of the radius of the core to the mass of the core, Ri/Mi, is such that Ri/Mi < 2.6 and that some of these 
modes can be excited by the perturbing spherical shell. If the phase transition is slow all of these models are stable (at 
least to radial oscillations). If the transition is fast, the star can oscillate radially. We derive an analytic expression 
for the square, w^, of the radial frequency of a relativistic DDS (see appendix A). We find that < for most (but 
not all) of the DDS stars whose w-modes are excited. The timescale of this instability is in many cases comparable 
to the period of the excited w-modes. While homogeneous stars have trapped modes only for 2.25 < R/M < 2.6, and 
the number of these modes decreases as R/M goes from 2.25 to 2.6, we will show that DDS models with R/M < 2.6 
and Ri/Mi > R/M can have, for some values of the jump 7, many more trapped modes than the homogeneous star 
with the same value of R/AI as the DDS. 

The paper is organized as follows: The metric coefficients and pressure of a static DDS star and a discussion of the 
way we parameterize these stars are given in Section ||. We then review, in Section [II , the perturbation model, a 
spherical shell surrounding the star, and we present the odd parity equation of motion of the matter in the shell, as 
well as the specific choice of multipole perturbation we use. In Section [V we present the wave equation governing 
odd parity perturbations due to the shell and review briefly the method we use to solve this equation. We then 
investigate in Section ^ the existence of trapped w-modes in DDS stars. In Section VL excitation of w-modes by 



the perturbing shell is cons idered and the stability of DDS is investigated in Section VII. Conclusions are briefly 
summarized in Section |VIIl| . The frequency of radial pulsation of a relativistic two-phase star with layers of constant 
and incompressible density is given in the appendix. 

Throughout the paper we use geometric units G = c = 1 (except where noted) and the conventions of Misner, 
Thorne and Wheeler Isl . 



II. THE BACKGROUND SPACETIME: DOUBLE DENSITY HOMOGENEOUS STAR 
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A. The metric coefficients 



We consider the equilibrium star to be static and spherically symmetric with a radius R and a mass M. The star 
has a inner core of radius r = Ri, mass Mi and constant density p+, surrounded by an outer layer of density p_ < p+- 
We denote the ratio of the two densities, by 7 = p+/p- > 1. When 7 = 1, the outer layer is absent and the star has 
constant density throughout its extension. The metric of this background spacetime can be written in the form 



sm 



Outside the star, the metric coefficients are 



-A(r) ^ I 



2M 



, r>R. 



(1) 



(2) 



Inside the star, we need to solve the hydrostatic equilibrium equations of general relativity (see for example 
equation (3)) under the assumption of two layers of constant density. This has been done by several authors (see e.g., 
pO,B5[). In this work we will use the metric coefficients in the form derived in p5|. 



and 
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m(r) = 



r<Rr 



(3) 
(4) 
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, R^<r<R 



R 



i? = l(l-37'^) (^l-47rp_^l-^ re3^W/2drj+^7-yi 
where the masses M^, M are (from ^) 

It is also useful to write down the expression for the pressure inside the star. 



2M 



P- 
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2M, 
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(5) 

(6) 
(7) 

(8) 
(9) 



(10) 



B. Specification of a double-density stellar model: Constraints on the values of R/M and Ri/Mi 

We will specify a stellar model by choosing the ratio of its densities, 7, the compactness of the star, R/M , and 
the compactness of its core, Ri/Mi. It is clear from the expressions (^,^ for the metric coefficients and ( [lO| ) for 
the pressure, that R/M and Ri/Mi must both be greater than 2. However this is not the only constraint on these 
parameters. Two additional constraints are that 7, R/M, Ri/Mi must be such that Ri/M exists and that the pressure 
(O) be positive everywhere inside the star. 
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Lets first address the constraints imposed by the condition that Ri/M should exist. To an arbitrary choice of the 
parameters 'y, R/M, Ri/Mi, corresponds three possibilities: a) There is no value of Ri/M; b) There is one value of 
Ri/M; c) There are two values of Ri/M. To understand this point, we notice that once "f, R/M, Ri/Mi have been 
specified, Ri/M can be determined from (||), 

'R,\' R, 7 [R^ V . 1 f RX" Q 



M J M^j~l\Mj 7 - 1 VM 

This cubic equation will have in general three distinct roots. From the theory of cubic equations we can show however 
that, if 

R MiV 4 7^ 



> i (-[2] 

M Rj 27(7-1)2 ' ^ ' 

two of the roots are complex and the third is negative. The only relevant set of parameters, 7, R/M, Ri/Mi, is then 
the one for which the inequality ( p^ ) does not hold. For that set, there will be three real roots. One of them is always 
negative and can thus be discarded. The other two, although positive, will only be mathematically admissible if 
Ri/M < Ri/Mi (since M > Mi). From an analysis of the two positive solutions we arrive at the following conclusions, 

• If 7 < 3, there is always a single solution Ri/M, if 

This condition states that the core should be less or as compact as the star. There is no solution otherwise. 
For given 7,i?/M, we have that Ri/Mi = R/M, corresponds to i? = Ri,M = Mi, i.e. the stellar model has a 
single phase with homogeneous density p+. As we then increase Ri/AIi above R/M we get ever smaller values 
of Ri/M which, in the limit Ri/Mi 00, reaches zero. 



If 7 > 3, there is again a single solution Ri/M, if 



Ri R /W.N 
M.'-M- ^''^ 

Otherwise, i.e. if the core is more or as compact as the star, 

Ri R 

M.^Tr ^''^ 

and if (^2|) does not hold, there will always be two solutions Ri/M. In particular if Ri/Mi = R/M, one of the 
solutions is the single phase star of density p+ as before (for which R = Ri,M = Mi), but a second solution 
generally exists with two phases and for which necessarily R ^ Ri, M ^ Mi. A point worth of mention is that 
for fixed 7, R/M, one of the two solutions Ri/M approaches zero as 7 ^ 00, while the other approaches Ri/Mi 
in the same limit. This is visible in figure 1, for the special case R/M = 5, Ri/Mi — 2.5. 



All the above considerations were derived solely from an analysis of the equation (|ll| ) . There might be cases however 
for which a stellar model is apparently possible from the above considerations but it is actually physically impossible 
because its pressure (^0|) is negative in some region of the star js^ . A necessary condition for positivity of the pressure 
everywhere inside the star, can be obtained by requiring that the pressure at r = Ri and at r = 0, be positive. This 
leads to the inequality, 



which is equivalent to 



1 _ 47rp_ ^1-^ y re3^('-)/2rfr > . (17) 

We can use either one of these inequalities to verify whether a set of apparently admissible parameters 
-f, R/M, Ri/Mi, Ri/M, characterizes a stellar model with positive pressure everywhere inside the star. If it does 
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not, we discard that set of parameters. Even if these inequahties are vahd however, there might be a small set of 
models for which the pressure is negative. Thus the best way to check whether a stellar model is possible is simply 
to plot the pressure at all points and see whether or not it is always positive. 

Figure 1 summarizes the results of this section. It shows the physical solutions Ri/M of (|ll]) as a function of 7, for 
fixed values of R/M, Ri/Mi. The lower curve, corresponding to the case R/M — 2.5, Ri/Mi — 5, shows that only one 
root exists for any 7 > 1. As 7 — > 00, the root approaches zero, i.e. the larger the jump in density the smaller the core 
is. If we switch the value of R/M with the value of Ri/Mi, making R/M > Ri/Mi, two roots are possible whenever 
7 > 51.94072268 (none exists otherwise). The two roots are equal to Ri/M = 1.6994 at this value of 7, but one of 
them approaches zero as 7 —> 00 (i.e. in this limit it approaches a single layer homogeneous star with density p-), 
while the other root approaches Ri /Mi (i.e. it approaches a single layer homogeneous star with density These 
limits can easily be derived from (pT|): in the limit 7 — *■ cx3, it reduces to a quadratic equation, with the solutions 
and Ri/Mi. The behavior of Ri/M as a function of 7 for other given values of R/M, Ri/Mi is qualitatively the same 
as the behavior displayed in figure 1. 

We would like to stress the following general property of DDS stars with R/M > Ri/Mi. Physical models of 
these stars only exist for 7 > jmin- The value jmin depends on the relative value of R/M and Ri/Mi. The bigger 
the difference between these two parameters, the bigger 7 is. For example for R/M — 5 and Ri/Mi = 2.5 it is 
7mm ~ 51.94072268, while for Ri/Mi = 2.3 and R/M = 5 it is 7™™ « 67.3019402, and for Ri/Mi = 2.3 and 
R/M = 6 it is 7™„ = 117.80657. 



III. MODEL: PERTURBATION BY MATTER MOVING ON A SPHERICAL SHELL 



The source driving the star away from equilibrium is a spherical thin shell of radius r — Rsheii > R, surrounding 
the star. Most of the details of this model have been presented in [|o) for the even parity case. We give in this section 
only those details that are relevant to odd parity perturbations. 

The stress-energy tensor of the shell is 



ST 



af3 
shell 



= Vl - 2M/rS"^S[r - Rsheii] • 



(18) 



We choose the matter in the shell, so that the only non-zero components of the surface stress-energy tensor, S*"^, 
have odd parity. The metric of the spacetime under this odd parity perturbation can then be written, to first order 
in the perturbation, as 



(0) , 



(19) 



where the "(0)" index denotes the background solution and e is the perturbation parameter. Because of the spherical 
symmetry of the shell, we can decompose the components SQ^,Sgip, of the surface stress-energy of the shell, in odd 
parity vector and tensor spherical harmonics. 



de 



Se 



fYio 
d9^ 



de 



(20a) 
(20b) 



Here Y/o is the scalar spherical harmonic. Upon substitution of (poj ) in the equations of motion of the shell (equations 
(8) of we get the single equation. 



hell 



dSlo 



(/-l)(? + 2) , 
1-2M/R,heii dt ^ 2 
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(21) 



By choosing 31^(1) we uniquely determine Sfg. In this work, we make the choice 



M 



(22) 



Since all perturbation equations will be proportional to e we will omit it henceforth. 
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IV. THE WAVE EQUATION GOVERNING ODD PARITY PERTURBATIONS 



For odd parity perturbations, all physical information can be obtained from a single function QiQ(r,t), that obeys 
the linear wave equation (see [gj for details), 

-Qp + U{1 + 1) — + 47rr [p{r) - p{r)] ^ — e ^ > = Sio{r, t) , (23) 

where r* is the usual tortoise coordinate, 

dr* ' ^ ' 

and the metric coefficients, pressure and remaining quantities characterizing the equilibrium DDS star where given in 
Section 0. The source term, constructed from the stress energy of the shell, has the form 



2D 



'5.0(0 = ^(1- — , (25) 



1 _ ^] (i _ ^] _ (i _ 



2 



2M 

Dio = SvrWf - Sfoimr - Rsheii] ■ (26) 

V ^shell 

We solve equation (^), by writing all time dependent quantities as Fourier integrals and imposing, at radial infinity, 
the condition that the wave be purely outgoing. At r* oo, we can then write, 

Qio{u = t-r*) = ^ f Aio{cu)e-"^"duj , (27) 

where 



M^Wl{cu)[l{l + f) - 2] V ~ Rshell 



y' {Rshell 7^) + Rshell ; {Rshell , ^) 

dr 



(28) 



Here y^'^^ is the regular solution of equation ( p3[ ) when the source term is zero and Wi{uj) is the Wronskian of the 
regular and outgoing solutions of (p3|). Details are given in We would like to point out that despite the fact that 
the density is discontinuous both at the inner core radius Ri and at the surface r = R, Qiq and dQio/dr should be 
continuous at both of these points. This follows from the requirement that the intrinsic and extrinsic geometries of 
the (odd parity) perturbed surfaces of the spacetime be continuous everywhere jS^] . 

In the remaining sections dealing with odd parity quasinormal modes, we will restrict attention to quadrupole 
{£ = 2) w-modes. 



V. THE TRAPPED ODD PARITY MODES OF A DOUBLE DENSITY STAR 

As is well known from the work of Chandrasekhar and Ferrari ||^, a homogeneous star has trapped odd w- modes 
whenever R/M < 2.6. The existence of these modes depends first on whether the potential, 

F.(.) = e^«{i(i±ll-^+4.(p-p)} , (29) 

of equation (|2^), has a well, and second on how broad and deep the well is. Homogeneous stars with R > 2.6M do 
not have a well and as such they have no trapped w-modes. 

In this section, we investigate the effect, on the trapped modes, of adding a second layer of constant density to an 
already homogeneous star. We have two reasons for undertaking this study. On one hand we want to confirm the 
argument given in Section |, that DDS stars with R/M > 5, do have trapped w-modes. As we show in the next 
section, some of these modes are excited by the perturbing shell. On the other hand we want to show that a DDS can 
have a very different w-mode spectrum from a single density layer stax, depending on the values of R/M and Ri/Mi. 

To compute the quasi-normal modes we have to solve the equation (p3) with a zero source term, a complex frequency 
and the boundary conditions of regularity at the center of the star and a purely outgoing wave at radial infinity. Only 
for a discrete set of complex frequencies (the odd parity quasi-normal modes of the star) will a solution to this 
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boundary value problem exist. The task of finding these frequencies is usually a complicated numerical problem. 
However, since we are only interested in the trapped w-modes and these are slowly damped, we can use the simpler 
resonance method [|8|-|o| to compute the frequencies. In this method, we integrate the homogeneous equation (|2^), 
for a given real frequency, starting at the center of the star and up to a point far away from the star. We know that 
the general solution of ( p3| ) at that far away point, is a linear superposition of an ingoing and an outgoing wave. At 
a quasi-normal frequency, the coefficient of the ingoing solution should drop abruptly to zero. By plotting say, the 
square of this coefficient, |yli„gp, versus the frequency, we can obtain, from the location of the minima, the real part, 
LOji, of the frequency of the least damped modes. Since near the minima, | Amgp = const[{ijj — lor)^ we can then 

compute its imaginary part w/. The method will only work if the imaginary part of the frequency is much smaller 
than the real part. 

Numerical experimentation suggests that a DDS star will have no trapped modes if 1) R/M > Ri/Ali > 3 or if 2) 
Ri/Mi > R/M > 3. And the reason is exactly the same as for a single layer homogeneous star: The potential ( [29| ) 
has no well in this case. We then restrict attention, in this section, to DDS models for which either 1) or 2) is not 
valid. Clearly an exhaustive presentation of results, even in the chosen range of R/M, Ri/Mi values is impossible. 
Instead we will present results for DDS with specific values of R/M, Ri/Mi and 7, that we believe illustrate most 
of the qualitative differences, as well as some similarities, between the trapped w-mode spectrum of a DDS and the 
spectrum of a single density star. 



A. Trapped modes of DDS stars with R/M > 5 

We turn our attention to the case that is one of the main motivations of this work: Do DDS stars with R/M > 5 
have trapped w-modes? A numerical investigation quickly shows that they will but only if Ri/Mi < 2.6. This means 
that all cases of interest for which R/M > 5, must have a core more compact than the star, i.e. Ri/Mi < R/M. As 
we recall from Section || (and figure 1), there will then be two Ri/M solutions of (11) as long as 7 > 7mm. None 



exists otherwise. These two solutions are equal for 7mm, but as 7 ^ 00, one, that we will call solution A, approaches 
Ri/Mi and the other, solution B, approaches 0. As we noted in the final paragraph of Section 7mm will depend 
on the relative value of R/M and Ri/Mi. It is higher the bigger is the difference between R/M and Ri/Mi. An 
important conclusion we can draw from all these considerations is that DDS stars with R/M > 5 and Ri/Mi < 2.6, 
require values of j > 46. 

To illustrate the existence of trapped w-modes for stars with R/M > 5, we compare in figure 2, three stellar models: 
A single density model with R*/M* = 2.3 and the two DDS models with R/M = 5, Ri/Mi = 2.3, 7 = 68, the one we 
called solution A, with Ri/M = 1.6442, and the one we called solution B with Ri/M = 1.4647. The choice 7 = 68 is 
very close to the minimum value of 7 (67.3019402) for which DDS models with R/M = 5, Ri/Mi = 2.3 exist and as 
such the two roots A and B are very close to each other. In figure 2a, we compare the potential (^9|) for each of these 
stellar models. The DDS potentials are nearly equal because the two roots Ri/M are very close to each other. We 
thus expect the DDS models to have a similar trapped w-mode spectrum, which we confirmed numerically. The DDS 
potentials are also deeper than the potential of the single layer star and hence we expect the double layer stars to have 
more trapped w-modes than the single layer star. That this is true is visible in figure 2b, where the least damped of 
these modes are shown both for the single layer star and for the DDS with Ri/M = 1.4647. While a single layer star 
has three or four trapped modes, the DDS star has seven or eight. The values of the frequencies of the two sets of 
w-modes do not appear to be easily related. For reference, we note that the three least damped w-modes of the single 
layer star have the frequencies w/M* = 0.235 + IQ-^z, wj/Af* = 0.316 + 5.6 x IQ-'^i and ujniM* = 0.392 + 7 x lO^^i. 
The six least damped frequencies of the DDS model are indicated in figure 2. The jumps that occur at r = Ri, R, in 
figure 2a, should be expected from ( p^ ) and the fact that the density is discontinuous at these points. It is also worth 
noting that while the maximum of the potential lies in the outer layer (of density p-), the minimum lies deep in the 
core (with density p+) of the DDS stars. 

What happens to the trapped w-mode spectrum of each of the DDS stars with R/M = 5, Ri/Mi = 2.3, when the 
jump 7 is increased and how does these spectra compare to the spectrum of the single layer star with R* /M* = 2.3? 
To understand what happens to the spectrum of solution A, as 7 ^ 00, we notice that this solution can be obtained 
from a star with R* /M* and density p+ by adding to it a layer of density p-, such as that the resulting DDS star 
has Ri/Mi — R*/M* — 2.3 and R/M — 5. As 7 ^ 00, Ri/M Ri/Mi (see figure 1), i.e. in this limit we again 
obtain the single layer star. Hence we expect that the w-mode spectra of solution A, will coincide with the w-mode 
spectra of the homogeneous single layer star for 7 — > cx). The two spectra should differ the most when 7 is near the 
minimum value allowed (in this case 67.3019402). The spectrum of w-modes of solution B behaves differently in the 
limit 7 ^ 00. This solution can be obtained from a star with R*/M* = 5 and density p^, by transforming a inner 
portion into a core of density in such a way that R/M = R* /M* — 5 and Ri/Mi = 2.3. As 7 ^ 00, Ri/M — > 0, 
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i.e. the inner core shrinks and we approach the original single layer star we started with. But a homogeneous star with 
R*/M* = 5 has no trapped w- modes (although it has w- modes). Therefore we should expect the trapped spectrum 
of solution B to disappear as 7 — s- oo. 

To investigate the behavior of the trapped w-modes for large 7, we compare in figure 3a, the potentials of the 
same three stellar models we used in figure 2, except that we consider a higher density jump for the DDS models, 
7 — 200. As we can see, the potential of the single layer star is very similar to the potential of the DDS solution 
A (with Ri/M — 2.179) and as such their least damped modes, shown in figure 3b for solution A, are very similar: 
ujjM = 0.221 + 4 X 10-6 (compare with 0.235 + lO'^), ujuM = 0.299 + 2 x W-^i (compare with 0.316 + 5.6 x 10""*^) 
and uJiiiM = 0.372 + 3 x 10"^^ (compare with 0.392 + 7 x lO^^i). On the contrary the potential of solution B (with 
Ri/M = 0.607) is very different from the other potentials and their least w-modes (shown in figure 3b) are therefore 
very different. Although all three models have four trapped w-modes with similar imaginary parts of their frequencies, 
the real part of the frequencies of solution B are significantly higher than the others as is clear from figure 3b. By 
increasing 7 well above 200, we confirmed that the solution B modes have real frequencies that approach infinite 
values, i.e. they disappear in the limit j ^ 00 and that the spectrum of solution A approaches the spectrum of the 
homogeneous star with R* /M* = 2.3 in the same limit, as we predicted above. A comparison of figures 2b and 3b, 
also shows that the A solution has more trapped modes the smaller 7 is (in 2b), and the reason why is evident from 
both figures 2a and 3a and the arguments above. 

From figures 2 and 3 we can conclude that DDS stars with R/M = 5,Ri/Mi = 2.3, and values of 7 near the 
minimum allowed, have more trapped modes than a single density star with R/M = 2.3. But as 7 is increased the 
number of trapped modes of the DDS stars quickly becomes the same as the number of single layer stars. Thus the 
effect of the second layer on the trapped modes seems to be more pronounced for the smallest possible values of 7. 

We only considered the case R/M = 5, Ri/Mi = 2.3, but DDS stellar models with Ri/Mi = 2.3 but higher values of 
R/M have trapped w-modes that behave similarly. The only main difference is that the higher R/M is (with Ri/Mi 
fixed), the higher the minimum value of 7 needs to be for the model to be physically possible, as we discussed in 
Section ||. This makes these DDS highly unrealistic as models of neutron stars with first order phase transitions. 

We also investigated models with Ri/Mi = 2.5, R/M = 5 and compared their trapped w-modes to the trapped 
modes of single layer stars with R/M = 2.5, which only have one such mode. Contrarily to what happens when 
Ri/Mi — 2.3, the DDS star never exhibits more than one trapped w-mode, no matter what the value of 7 is. The 
frequency of this w-mode however, is highly dependent on the value of 7. This dependence can be understood using 
the same reasoning we presented above. For Ri/Mi > 2.6 and R/AI > 5 the DDS stars seem not to have trapped 
modes. 



B. Trapped modes of DDS stars with R/M < Ri/Mi 

While DDS stars with R/M > 5 can have a very different w-modc spectra from homogeneous single layer stars, it 
is for DDS with R/M < Ri/Mi that the difference between the modes of a two layer star and a one layer star is more 
striking. 

Trapped modes exist only for R/M < 3, so we choose as illustrative examples the DDS stars with R/M — 
2.3, Ri/Mi = 5 and R/AI = Ri/Mi = 2.5. The interesting aspect of these stars is that, for relatively low values of 7 
(presumably as small as a few percent) the w-mode spectrum can be extremely different from the w-mode spectrum 
of a single phase star with the same value of R/M. We should keep in mind that, from Section |l| (and figure 1), if 
R/M < Ri/Mi, there is only one solution Ri/M. This solution approaches zero as 7 ^ 00. 

We thus consider a single layer star with R* /M* — 2.3 and density p_. The three least damped w-modes of this 
star, are visible in figure 2b and were discussed above. We then transform a inner layer of this star into a core of 
density p+ and Ri/Mi = 5, in such a way that R/M = R* /M* = 2.3. As we recall from Section there will be in 
general one solution Ri/M of ( |lT| ) for any value of 7. For R/M = 2.3, Ri/Mi — 5 however, models with 7 < 17.6 
have negative pressures near the core and are thus not allowed. For 7 > 17.6, the value of Ri/M is maximum for 
the smallest allowed value of 7, and it decreases to zero as 7 ^ 00 (see figure 1). In other words, the configuration 
R/M = 2.3, Ri/Mi = 5,7 00, is again the single layer star with density p_, that we started with. This means that 
the w-mode spectrum of a R/M = 2.3, Ri/Mi = 5 star, should differ the most from the spectrum of the R* /M* — 2.3 
star, for values of 7 near 17.6. The two spectra should approach each other as 7 — *■ 00. This is shown explicitly in 
figure 4. 



Figure 4a, shows the potential (29), for three stellar models: The single density star with R*/M* = 2.3 and two 
R/AI = 2.3, Ri/Mi = 5 DDS stars, one with 7 = 18, the other with 7 = 28. Figure 4b allows us to visualize the real 
part of the trapped modes of these models and thus compare the effect of the additional layer. Figure 4b, clearly 
shows significant differences in the three spectra of oscillation modes. These differences can be partially understood 
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by looking at the shape and depth of the potential well for each of the three models, in figure 4a. There, we see that 
the smaller the value of 7 is, the broader and deeper the potential well is. We thus expect the number of trapped 
w-modes for the model 7 = 18, to be higher than for the model with 7 = 28, and this in turn to be higher than the 
number of trapped modes of the single layer star. That is precisely what is observed in figure 4b. For 7 = 18, there 
is a very large number of trapped w-modes. For 7 28, however only a few (about 10) remain and only three or four 
for the single layer star (which is the limit 7 00). We should also note from figure 4a, that while the maximum of 
the potential always resides outside the star (in contrast with the stars in figures 2a and 2b where it resides in the 
outer layer), the minimum lies deep in the inner layer for the model 7 = 18, but lies at the boundary core-envelope 
for the model 7 = 28. 

Figure 4 focused on a particular case, for which R/M < Ri/Mi. It is clear however, that the reasoning we used 
can be applied, with qualitatively similar results, to stars with R/M — 2.3, but higher values of Ri/Mi. The only 
significant change is that as the core becomes less compact (i.e. as Ri/Mi increases), the minimum allowed value of 
7, for which a solution of (^l|) exists and has a positive pressure everywhere, decreases. While for Ri/Mi = 5 it was 
« 17.6, for Ri/Mi > 7 all values of 7 > 1 are allowed. 

Our final example, considers the DDS with R/M — Ri/Mi = 2.5. This case is interesting because, depending on the 
value of the density jump, 7, there can be two, not just one, DDS stars with different values of Ri/M. Thus, for any 
value of 7, one of the solutions of ( [ll| ) is Ri/M — 2.5, i.e. Ri — R, which is just a single layer star with homogeneous 
density p+ throughout. For 7 < 3 this is the only solution. For 7 > 3, however, there is a second solution with two 
layers and Ri < R. This solution has however negative pressure in the inner core, for 3.87015 < 7 < 48.374. For 
this range only the single layer star solution exists. But in the ranges 3 < 7 < 3.87015 and 7 > 48.374, the second 
solution has positive pressure everywhere and we do have two different Ri/M solutions with R/AI = Ri/Mi = 2.5. 
For 7 = 3, the two solutions are the same. In the limit 7 — > 00, the second solution corresponds to a single layer star 
with density p_ and R/M = 2.5 since Ri/M (see figure 1). It is then clear that the trapped w-mode spectrum 
of each star, in these two limits, will be the same and equal to the w-mode spectrum of a single density star with 
R/M = 2.5. But for values of 7 in between these two extrema, the stars modes are very different. 

In figure 5a we show the potential (|29|) , for stellar models with R/M = Ri/Mi = 2.5 but four different values of 7: 
3 (single layer star with Ri = R), 3.5 (a double layer star with Ri/M = 2.158), 49 (a DDS with Ri/M = 0.388) and 
200 (a DDS with a small core, Ri/M = 0.184). The models are identified in the figure by their Ri/M values. In figure 
5b we show the least damped w-modes for the four models depicted in figure 5a. The striking differences can again be 
understood from figure 5a and the discussion above. Thus the DDS model with 7 — 200 has a single trapped mode, 
whose value is nearly equal to the single trapped mode of the single layer star with 7 = 3. This is expected since the 
DDS is very close to be a single layer star. If we compare the DDS model with 7 = 3.5 to the single layer model with 
7 = 3, we notice that a second trapped w-mode appears. And many more appear for 7 = 49 (whose potential is the 
deepest in figure 5a). 

We have also studied DDS models with R/M = Ri/Mi = 2.7 and R/M = Ri/Mi = 2.8. These are similar to the 
case R/M — Ri/Mi — 2.5, in the sense that there is, besides the single layer solution, a second solution with two 
layers. We found that for certain values of 7, this second solution has one trapped mode. This should be contrasted 
with the single density homogeneous solutions with R/M — 2.7,2.8, which have no trapped modes. Indeed as was 
shown in these stars have trapped modes only for R/M < 2.6. 

One of the conclusions of this section is that the w-mode spectrum of a DDS star can differ greatly from the 
spectrum of a single density star. The differences are more pronounced in the (physically less realistic) case in which 
R/M < Ri/Mi. Indeed in this case, we can form DDS models with a much higher number of trapped w-modes than 
the single layer star trapped modes as is shown in figures 4b and 5b. 

The other conclusion, is that DDS models with R/M > 5, have trapped modes as long as Ri/Mi < 2.6. This 
requires 7 > 46, which might be hard to attain in the universe. It is interesting to note that the more significant 
difference between double and single density star's trapped mode spectra happens, for values of 7 near the minimum 
allowed. As we will show in the next section, some of the w-modes of these stars can be excited. 

VI. EXCITATION OF W -MODES OF STARS WITH R/M > 5 

To illustrate w-mode excitation, by matter moving in the spherical shell surrounding the star, we will consider three 
DDS models and one single density model. We choose the models so that all have the same value of core compactness, 
Ri/Mi = 2.5. We also choose in all cases, the radius of the shell to be Rsheii = 50 and the Gaussian parameter in 
( p2| ) to be aM^ = 0.1. Figures 6 and 7 show the waveform for each of these four models. The general behavior of the 
waveform is similar in all four cases: Since the shell is located at r = 50, in all waveforms there is a sharp pulse at 
u ~ — r*(50) « —56 corresponding to the motion of the matter in the shell at t « 0. Part of this pulse travels inward 
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and is reflected by the star. The reflected pulse, whose form and shape depends on the star, will appear at u « 56. 
There is then a smafl transient period, followed usually by ringing at one (or more) of the w-mode frequencies of the 
star. 

In figure 6a we present the waveform for a star with homogeneous density p+ throughout and R/M — 2.5. The 
least damped (odd parity) mode of this star is loM — 0.412 + 0.022i (barely visible in figure 5b) and is clearly excited 
after u « 93. 

In figure 6b we consider the DDS star with R/M — Ri/Mi = 2.5,7 — 3.5. As we discussed in there are 
two solutions for these parameters. One is the single density star considered in figure 6a. The other is a DDS 
star with Ri/M = 2.15831. This star has two clearly visible trapped ui-modes in figure 5b. Their frequencies are 
LOiM = 0.268 + lO^^i and lohM = 0.357 + 0.003i. The third w-mode (barely visible in figure 5b) has frequency 
ujjjjM — 0.412 + 0.05i. In figure 6b the star starts ringing at the frequency ivjjjM at m « 130. At m « 170 there is 
an interference between modes w/// and w// and after u « 270 the star rings at the frequency cj//. The least damped 
mode does not appear to be excited. 

In figure 7 we consider stars with R/M > 5. We also choose 7 to be very close to the minimum value allowed for the 
values of R/M, Ri/Mi indicated. For R/M = 5, Ri/Mi = 2.5,7 = 52, there arc two possible models: Ri/M = 1.6666 
and Ri/M = 1.732. In figure 7a we consider the first of these models. Ringing of the least damped mode, with 
frequency ujM = 0.482 + i0.014, sets in after u w 100. In figure 7b it is the DDS model with R/M = 6,R^/AI, = 
2.5, Ri/M = 1.598,7 — 92 that is considered. Ringing at the least damped frequency {ujM — 0.533 + iO. 021) sets in 
after u « 100. 

Two points are worth noting: 1) As is clear from figure 6, the gravitational wave emitted by a single layer star 
might be considerably different from the wave emitted by a DDS with the same degree of compactness (in this case 
both stars in figure 6a and in figure 6b have the same value of R/M), since the two stars have very different slowly 
damped quasi-normal modes. The interesting point is that this is true even for small values of the density jump, 
namely 7 = 3.5. But only for very compact stars (i.e. when both R/M and Ri/Mi are < 3) will this likely occur. 2) 
DDS models with R/M > 5 can have their w-modes excited. But this happens at a price: The core Ri/Mi < 2.6, 
which might be a somewhat realistic assumption for neutron stars, but the jump 7 > 46, which seems to be unlikely. 



VII. INSTABILITY OF DDS STARS WITH A FAST PHASE TRANSITION 



While DDS stars with R/M > 5 might have their least damped u;-mode(s) excited, it is important to ask whether 
these models are stable over timescales greater than the period of the excited z«-mode(s). To address this question 
it is useful to divide the DDS stars into two extreme cases: DDS with a very slow phase transition and DDS with a 
very fast phase transition. We model the first case as if no phase transition occurs and we model the second case as 
if the phase transition is instantaneous. Although we might question the validity of this last assumption, we recall 
that the constant density assumption behind the DDS models is also unrealistic since it allows infinite sound speeds. 
Nevertheless we can obtain useful limits on the stability behavior of stellar models with a phase transition and it is 
with this in mind that we use the two limits of slow and fast phase transitions. 

The adiabatic index on either side of the interface between the two phases is infinite since the phases are supposed 
to be incompressible and of homogeneous density. In more realistic models of stars with phase transitions, the two 
phases will have different and finite adiabatic indices. For those models it might be possible to obtain stability criteria 
in terms of the values of the adiabatic indices, in analogy with what is done in single phase stars. 

No matter what the speed of the phase transition, a nonrotating (no £ — 1 perturbation allowed) DDS star is 
stable to odd parity perturbations. This can be easily shown using the fact that all these perturbations (with £ > 2 
multipoles)must obey the wave equation (p3|). From (p3|), an "energy-like" conservation law for the perturbation can 
be obtained which implies that the perturbation cannot grow without bound at any r, i.e. that the star is stable to 
odd perturbations. For the original proof see Also |Q. 

Stars with a very slow phase transition (slow compared to the periods of the relevant w-modes) are dynamically 
stable to radial and probably any even parity perturbation. But stars with a fast phase transition, might be unstable 
to even parity radial oscillations as shown in |2^J23] (Newtonian models) and in 



24 , 2a] (general relativistic models) . 



Indeed in appendix A, we derive the frequency of such oscillations, and from it we conclude that the star is unstable 
whenever < 0, i.e. whenever. 
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a result obtained by Ovakimova ||2^ . 

For very small cores, (i.e. i] or Ri << R) the term in square brackets is dominant over the other two terms, 
and we conclude that the DDS star will be unstable whenever its density jump is 

3y/T^2M/R_ 1 

7 > t: — = - — jj , (31) 

2 - 2M,/i?, 1 - Anp^^l^2M/Rjl^^ re^^/^dr 

as was noted first by Seidov |2^. In the Newtonian limit (R/M, Ri/Mi << l,p_A/^ << 1) this reduces to the well 
known result ||2l|,^ 7 > 3/2. 

We will now use (AI6) to compute the frequency of radial oscillation and to find out when are the stars with a 
fast phase transition unstable. We chose to display results (in figure 8) as plots of the frequency sgn{uj'^)\ijjM\ (not 
uj^M"^), where sgn{uj'^) — ± depending on the sign of o;^, versus the density jump, 7. For each plot we keep R/M and 
Ri/Mi constant. In this way, when the star is unstable we have negative numbers in the y axis. The inverse of these 
numbers multiplied by 2n gives the instability time of the star. The purpose of figure 8 is to give a rough idea of 
when stars with a fast phase transition are stable and when are they unstable. Even if the model is unstable however, 
it might be so on a timescale much longer than any of the w-modes that are likely to be excited. Therefore we will 
compare, for a few stellar models, the instability time with the damping time of the least damped i(;-mode(s) to get 
an idea of the importance, or not, of this radial instability. 



A. The case R/M > Ri/Mi 

In this case there are in general two solutions of ( pT| ) (see figure l)for 7 > jmim one with a small core and the other 
with a big core. We represent in figure 8a, three families of models with fixed values of R/M, Ri/Mi but different 
values of 7. The small core solution is represented by a dashed line, and the big core solution by a solid line. The 
frequencies are monotonous functions of 7 for each of the solutions. We see that the small core solution is always 
unstable for models with Ri/Mi < 2.6. The solution with the big core might be stable if 7 and Ri/Mi are high enough, 
(for example almost all models with Ri/Mi — 2.3 are unstable, but many of the models in the second solution with 
Ri/Mi = 2.5 and high enough 7 are stable). In particular all models considered in Section ^ arc unstable. Is the 
instability relevant in the context of w-mode excitation? Let us consider the model R/M = 5, Ri/Mi = 2.5,7 = 52 
and Ri/M = 1.67 (figure 7a above). The excited w mode has a damping time of tdamp/M = 27r/0.014 « 449. The 
instability time is 27r/0.072 « 87.3. For a star with one solar mass, these times are idamp — 2.2msec, tinst — OAmsec. 
Therefore for DDS stars with these parameters, and a fast phase transition, it is not relevant to talk of w-mode 
excitation. A similar computation for the model in figure 7b leads to the same conclusion: The damping time of 
the excited w-mode is tdamp/M « 299 and the instability time is tinst/M = 27r/0.06 ~ 105. For a star with one 
solar mass, these times are tdamp ~ ^■5msec,tinst = 0.5msec. All models with two layers in figures 2 and 3 are also 
unstable. 

But there are stable models: For R/M = h, Ri/Mi — 2.6,7 = 77, the small core solution {Ri/M = 1.005) is 
unstable, but the big core solution {Ri/M = 2.331) is stable. In general however, all solutions with R/M > 5, 
Ri/Mi < 2.6, for which the likelihood of u;-mode excitation seems the highest, will be unstable for the least possible 
values of 7. 



B. The case R/M < Ri/Mi 

As we remember from Section |l|, if R/M < Ri/Mi, there is only one solution of (|ll|). In figure 8b we show some 
curves frequency versus 7 for a few stellar models with quite small value of Ri/Mi. Th e for m of these curves is similar 



for all stars with R/M < Ri/Mf. for small 7 (near 1), the term in square brackets in (A16) completely dominates the 
other terms and the frequency has its maximum value. This value is positive. As 7 increases, the frequency drops, 
becoming negative for some models or staying positive for others, until eventually it reaches a minimum, after which 
it rises again approaching an asymptotic limit as 7 00. In figure 8b we consider three curves corresponding to 
DDS models with the same value of Ri/Mi = 5 but different values of R/M: 2.3, 2.5 and 4. We also consider a curve 
corresponding to models with Ri/Mi = 2.3 and R/M — 50. The three curves with Ri/Mi = 5, show that above a 
certain value of the jump 7, all the models are unstable. The frequency decreases as the jump in the density increases 
meaning that the instability time of these models increases. For small values of 7 the instability time is a minimum 
though. And this minimum can be much smaller (by several orders of magnitude) or comparable to the damping time 
of the w-modes. The conclusion is that models with R/M < Ri/Mi, having trapped modes will be certainly unstable 
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for almost all small values of 7 and that instability certainly kills the star before significant w-mode ringing goes on 
for small values of 7. However for very high values of 7, the instability might allow w-mode ringing. 

It is instructive to look at some of the models we considered in Sections |^ and Vl . The model considered in figure 
6a is obviously stable (it is a single density homogeneous star). The model in figure 6b is also stable. The remaining 
models, all with R/M — Ri/Mi — 2.5, in figure 5 are both unstable. All the models in figure 4 (with the exception 
of the homogeneous model) are unstable, as is evident from figure 8b. 



VIII. CONCLUSIONS 



We have shown in Section ^ that adding a second layer to a homogeneous density star, can change completely 
the w-mode spectrum. In particular stars with R/M > 5, very compact cores, Ri/Mi < 2.6 and high density jumps 
7 > 46, can have trapped (very slowly damped ) w-modes. Some of these modes are excited as we showed in Section 



VI. But where the effect of a second layer is more striking is when both R/M and Ri/Mi are smaller than 3. In that 
case, even for very small values of the density jump (7 ~ 18 in figure 4 and 7 ~ 3 in figure 5) we might have stars 
with completely different w-mode spectra from the spectra of homogeneous stars. We have found that DDS models 
with R/M = 2.5, can have a very large number of trapped w-modes. By contrast, homogeneous stars with a very 
large number of trapped w-modes, require values of R/M w 2.25, very near the maximum compactness of the star. 

Most of the DDS models with trapped modes, the most likely to have their w-modes excited, are unstable to radial 
oscillations, if the phase transition, at the origin of the density jump, is very fast. But if the transition of one phase 
to the other is slow, the models are stable. It would be of interest to investigate whether, in realistic stars, the phase 
transition happens on a slow or fast timescale, since this might determine whether such stars can or can not exist. 

Although we showed that DDS stars with R/M > 5 can have their w-modes excited, all of these stars require 
relatively large discontinuities of the density. Values that are unlikely to occur in realistic neutron stars. It is thus 
worthwhile to use different, more realistic, equations of state with a density discontinuity and see if the phenom- 
ena observed in this work for relatively high 7, occurs for realistic values of 7. For example does dropping the 
incompressibility and/or the constant density hypothesis change in a significant way w-mode excitation? 

It is also useful to ask whether the results reported here depend on the existence of a density discontinuity. We 
know that in a realistic neutron star, the density is usually continuous but might change over very small length scales. 
Thus it would be interesting to consider what would happen if the discontinuity is smoothed out and instead we 
consider the jump from p+ to p_ to happen continuously but over a very short length. Would it be possible, to have 
trapped w-modes and significant w-mode excitation, for stars with continuous density throughout and R/M > 5? 
An obvious idea based on this work, is to consider stars with two regions of very different densities (a inner higher 
density layer, surrounded by an outer low density layer), connected by a thin shell where the density drops abruptly, 
but continuously, from its high to its low value. Will the w-mode spectrum and excitation depend on the detailed 
form of the transition from high to low density? If yes, what are the relevant factors? 

We restricted attention in this work, to odd parity perturbations. It is of some interest to study even parity 
perturbations of DDS stars and see whether the excitation of w-modes of stars with R/M > 5 and Ri/Mi < 2.6 is 
comparable to the fluid mode excitation. We have started addressing some of these questions, in the hope they might 
shed some light on the issue of the relevance of w-mode excitation in realistic stars. 
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APPENDIX A: FREQUENCY OF RADIAL OSCILLATION OF A RELATIVISTIC TWO-PHASE STAR 

We derive the frequency of radial oscillation of a star with two-phases, each of constant density and under the 
assumption that the fluid is incompressible in each phase. We also assume that the phase transition is so fast that at 
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any instant during the radial vibration, the pressure on either side of the phase boundary is equal to the pressure at 
which the two phases are in equilibrium. 

Under these assumptions, Grinfeld and Bisnovatyi-Kogan and Seidov |^^, obtained the frequency of radial 
oscillation of a two-phase Newtonian star and from it the stability criteria previously derived by Ramsey . Haensel 
et al 1^ , Bisnovatyi-Kogan and Seidov |Q and Migdal |Q have also computed the radial frequency of oscillation of 
Newtonian two-phase stars under more realistic assumptions. The frequency of a relativistic two-phase star appears 
not to have been derived before, although Ovakimova obtained the criteria of instability of these stars using a 
method that did not require the knowledge of the frequency. 

In the following, let ^ denote the radial displacement of a fluid element from its equilibrium position and let the 
quantity 6h, be defined by 

Sh = (Al) 
P + P 

where Sp is the Eulerian variation of the pressure, and p, p are the equilibrium values of the pressure and density. 
A superscript "— " means the quantity is evaluated aX r = Ri — e and a superscript means it is evaluated at 
r = Ri + e {e « 1). It is straightforward to show, under the above assumptions, that the quantities 5h and ^ must 
obey the following three boundary conditions: 

1. Continuity of 5h at the unperturbed phase boundary 

6h+{Ri)^dh-{Ri). (A2) 

2. Continuity of the Lagrangian perturbation of the pressure at r = Ri, 

[p{R.) + p+]6h- +r^ = [P{R^) + P-]5h+ + t%r- (A3) 



Using (A2) and the equation of hydrostatic equilibrium to eliminate dp^ /dr, this equation reduces to 

VI, + 4T:Rfp{l 
R,{R,^2M, 



(p+ - p^)5h = {r{R^)[p+ +P{R^)] t{R.)[P- +P{R^)]] ^itp.^^^ll^;^ ■ (A4) 



3. Vanishing of the Lagrangian perturbation at the surface of the star 

We should note that ^+(i?i) ^ £^~{Ri), because there is a flux of mass across the phase boundary due to the phase 
transition. This should be contrasted with other studies [ ^3p^ ] of oscillations of stars with density discontinuities but 
without phase transitions, for which ^~^{Ri) = (Rj). 

Imposing the three boundary conditions ( A^, A4 , A^) to the solutions of the equations of radial oscillations, leads 



directly to an analytic expression for the frequency of such oscillations, as we now show. We will follow the equations 
and notation in chapter 26 of p^ ] in the remainder of this appendix, taking care of the difference between the form 
of the metric (|^) we use here, and the form they use. 

The Eulerian perturbation of the density is given by equation (26.11) of |3^. Noting that the star is incompressible, 
this equation reduces to 



J.2 fj^j. 



5X 
2 



= 0. (A6) 



But from equation (26.15) of |Q, we know that 

SX = ^{Vr + X,r)C- (A7) 

Substituting this equation in (|A6|), we get the following equation for ^, 



^__d 

^Jr'' ^' 2 



whose integral is (.{r,t) — A{t)e'^'^^^^'^ /r^ . Since ^ must be regular at the origin we have 
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r < Ri 



A{tY^^, R,<r<R. 



(A9) 



To determine 5h, we need the relativistic Euler equation, ( |Q, equation (26.18)), 

1 



and 



5v r 



6X + 87rr(p + p)6he^ 



(AlO) 



(All) 



which is obtained by c ombi ning equations (26.13b ) and (26.14) of B 



Substituting Svr in (AlC), by its expression in (|All|), (and using also A7) we obtain the following equation for Sh, 



1 



- 47rr(p + p)6he^ 



(A12) 



Let us assume that all quantities have a time dependence of the form e"^* where uj is the frequency of radial oscillation. 
Upon substitution of the expression ( [A9| ) in (A12), we obtain, after performing the radial integration, 



r<R^ 



Aj;^ dr {^2e3A/2j^ + + X,r)^e-^+^/^ [i.,. + ^]} + C\ R, < r < R. 



(A13) 



If we now impose the three boundary conditions (A2,A4 A5) we can determine to directly. The result is 



1 



'A 



i?3 i? - 2M 2 



Ri 



V.r + 



1 



B ^'^(R^)+iKR.)/2 AI, + 4TTRfp(RA 

T = [P[Rt)+P-\ 54 ■ 

A P+- P- Rj 



(A14) 



(A15) 



Using the metric coefficients in Section |l| and the general relativistic equations of hydrostatic equilibrium, we can, 
through lengthy but straightforward algebra, rewrite the frequency of radial oscillations in the somewhat simpler 
form, 



47r/9_ 



1 



1 



jl^^dre^^/^/r^ [R^^1-2M/R 



ry[(7 - 1)??' + 4] + 



1 - 2M/R 



i?,(7-l)(l-2Af,/i?,)3/2 
_36^p_(l-^) Tdre^^W/^ 

R J R 



y/1 - 2M/R V Rt JR 

where 77 — R J R is the ratio between the radius of the core and the radius of the star 



(A16) 
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FIG. 1. Roots Ri/M of equation (|ll|) when R/M = 5,Ri/A'U = 2.5 or vice-versa, R/M = 2.5, Ri/Mi = 5, plotted as a 
function of 7. 
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FIG. 2. Figure 2a displays the odd parity potential for 1 = 2 and three stellar models: a single density model with 
R/M = Ri/Mi = 2.3 and two double density models both with 7 = 68, R/M = 5, Ri/Mi = 2.3 but Ri/M = 1.4647 
and Ri/M — 1.6442. Figure 2b shows the least damped ?j;-modes of the single layer star and the double layer star 
with Ri/M = 1.4647. The six lowest frequencies of this double star are loiM = 0.180 + 10~*i, loijM = 0.245 + 10~^i, 
uJiiiM = 0.309 + 29 X 10"S, ujivM = 0.373 + 15 x 10"^j, ujvM = 0.434 + IQ-^i, uviM = 0.493 + 5 x 10"^j. 
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FIG. 3. Figure 3a displays the odd parity potential for ^ = 2 and three stellar models: a single density model with 
RjM = RijMi = 2.3 and two double density models both with 7 = 200, i?/M = b,Ri/Mi = 2.3 but Ri/M = 0.607 
and Ri/M = 2.179. Figure 3b shows the least damped lu-modes of the two DDS stars. The four least QN modes of the star 
Ri/M = 0.607 are ujiM = 0.524+3 x 10"^i, unM = 0.709+18x 10~^ lj/z/A/ = 0.886+3 x 10~^i and uivM = 1.054+2 x lO^^i. 
The frequencies of the other star are indicated in the text. 




FIG. 4. Figure 4a displays the 1 = 2 odd parity potential for three stellar models: A single density star of 
R/M = 2.3, a double density star with 7 = 18, ii/M = 2.3, /Mi = "a, Ri/M = 0.3817 and a double density star with 
7 = 2?,, R/M = 2.3, Ri/Mi = 5, Ri/M = 0.304. The real part of the least damped quasi-normal modes of these three stellar 
model are the minima in the graphics of figure 4b. The double density stars are identified, in figure 4a, by their Ri/M value 
and the single density star by its R/M value. 




FIG. 5. Figure 5a displays the 1 = 2 odd parity potential for four stellar models, all with the same R/M = Ri/Mi = 2.5, 
but difi^erent values of 7: A single density star with Ri/M = 2.5, and three double density stars with 7 = 3.5, Ri/M = 2.158, 
7 = 49, Ri/M = 0.388 and 7 = 200, Ri/M = 0.184 As noted in the text, for the three models with 7 > 3, there is always the 
solution corresponding to Ri/M = 2.5, in addition to the solution displayed here. Figure 5b shows the least damped lo-modes 
for each of the four models. 
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FIG. 6. Quadrupole {£ = 2) waveforms for two stars with R/M = 2.5 perturbed by matter in a shell located at Rsheii = 50M 
and with a Gaussian parameter aM^ = 0.1. In figure 6a, a homogeneous star is considered. In figure 6b, a DDS with 
Ri/Mi = 2.5,7 = 3.5 and Ri/M = 2.158. 
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FIG. 7. Quadrupole = 2) waveforms for two stars with Ri /Mi = 2.5, perturbed by matter in a shell located at Rsheii = 50M 
and with a Gaussian parameter aM^ — 0.1. In figure 7a, a DDS with R/M = 5,7 = 52, Ri/M = 1.67. In figure 7b, a DDS 
with R/M = 6, 7 = 92 and i?»/M = 1.598. 
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FIG. 8. Figure 8a displays the radial frequencies of oscillation of a two-phase star. The arrow labeled "line 1" represents 
the frequencies of DDS with R/M = 5, Ri/Mi = 2.5. The upper curve (in solid) is the solution that approaches Ri/M = 2.5. 
The lower curve (in dashed) is the solution Ri/M that approaches 0. The arrow labeled "line 2" represents the frequencies 
of DDS with R/M = 5, Ri/Mi = 2.3. The upper curve (in solid) is the solution that approaches Ri/M = 2.3. The lower 
curve (in dashed) is the solution Ri/M that approaches 0. The arrow labeled "line 3" represents the frequencies of DDS with 
R/M = 6, Ri/Mi — 2.3. The upper curve (in solid) is the solution that approaches Ri/M = 2.3. The lower curve (in dashed) 
is the solution Ri/M that approaches 0. Figure 8b considers models with fixed values of R/M and Ri/Mi > R/M, but varying 
jump 7 in the density. All the models are unstable for 7 > 5 but there is a small region near 7 1 of stable models. 
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